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Abstract 

Recent theoretical studies have predicted magnetic states for Ru and Os 
in the body centered tetragonal structure (bet) with c/a < 1. In this study, 
we present first principles calculations of the phonon dispersion for ferromag¬ 
netic Ru- and Os-bet along the epitaxial and uniaxial Bain paths, to eval¬ 
uate their dynamical stability. The phonon dispersions were computed us¬ 
ing the density functional perturbation theory, including the gradient correc¬ 
tions to the exchange-correlation functional within the plane-waves ultrasoft- 
pseudopotential approximation. The phonon dispersion for the local mini¬ 
mum in the Bain path with c/a < 1 as well as the uniaxial and epitaxial 
strained structures are analyzed. We find imaginary frequencies along differ¬ 
ent directions of the Brillouin zone, which indicates that both systems are 
dynamically unstable. Consequently, ferromagnetic Ru and Os in the body 
centered tetragonal structure with c/a < 1 are not truly metastable phases. 


1 Introduction 

The search for metals in metastable phases is a common approach to finding new 
materials. A metastable phase corresponds to a state in local mechanical equili¬ 
brium with a crystalline structure other than its ground-state, and it is therefore 
energetically less favorable [1]. Thus, one of the most active research helds in com¬ 
putational material science is the search for metastable structures with enhanced 
physical and chemical properties. Among these properties, ferromagnetism is one of 
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the most important properties of matter because of its technological relevance and 
the low number of ferromagnetic elemental metals. Because the crystalline structure 
plays a major role in the dehnition of the properties in any material, a possible way 
to induce ferromagnetism in metals with paramagnetic ground-state is by changing 
its crystalline structure to reach a metastable state. 

Ruthenium (Ru) and Osmium (Os) are transition metals with hexagonal close 
packed structure (hep) and paramagnetic ground-state. Nevertheless, Kobayashi 
et al. [2] studied bulk Ru in the body-centered cubic structure (bcc) by means of 
hrst principles calculations using the APW method and the LDA approximation 
for the exchange-correlation functional, and they hnd ferromagnetism for 5% of the 
lattice constant expansion. Motivated by this result, Siiki and Hio [3] experimen¬ 
tally studied the formation of Ru thin-hlms over substrates with different lattice 
constants, in an attempt to obtain ferromagnetic bcc Ru. Instead of obtaining the 
Ru-bcc structure, they hnd that the X-ray and electron beam diffraction patterns 
for Ru thin-hlms over Mo(llO) suggest the formation of a body centered tetragonal 
structure (bet) with lattice constants of a = 3.24 A and c = 2.69 A (c/a = 0.83). 
Subsequently, based on hrst principles calculations using the FLAPW method and 
the GGA approximation for the exchange-correlation functional, Watanabe et al. 
[1] studied bulk Ru in the bcc, hep, and face-centered cubic (fee) structures as well 
as its tetragonal deformation path, which is better known as the uniaxial or classical 
Bain path. They determined that the total energy in the Bain path revealed a local 
minimum in the bet structure for c/a = 0.83 and a maximum in the bcc structure, 
which indicated that the Ru-bcc structure is unstable. Recently, Schonecker et al. [5] 
using the FPLO method and the GGA approximation for the exchange-correlation 
functional, performed density functional theory (DFT) calculations of the epitax¬ 
ial Bain path for 28 non-magnetic ground-state transition metals to search for bulk 
ferromagnetic solutions assuming that thick dims could be coherently grown on sub¬ 
strates with fourfold symmetry. They find that Ru and Os in the bet structure with 
c/a < 1 presents a ferromagnetic order. Furthermore, Odkhuu et al. [6] employing 
DFT calculations with the FLAPW method and the GGA approximation for the 
exchange-correlation functional, predicted a large perpendicular magnetocrystalline 
anisotropy for bulk Ru in the bet structure, which is two orders of magnitude greater 
than conventional magnetic metals. These reports show the growing interest in the 
magnetic properties of Ru and Os in the bet structure. 

It is important to note that in bulk metals, a crystalline structure could corre¬ 
spond to a local minimum in the total energy along the Bain path, even though it 
could be dynamically unstable in that minimum [1]. Moreover, based on DFT cal¬ 
culations it has been reported that cubic monoatomic metals are elastically unstable 
in the bet minimum of the Bain path [7], which suggests that this result could also 
be the case for the hexagonal close packed metals. Regardless of the prediction of a 
ferromagnetic order in the local minimum of the Bain path for Ru and Os in the bet 
structure with c/a < 1, a study on the dynamical stability in these materials has not 
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yet been reported to the best of our knowledge. From a theoretical point of view, 
the stability analysis of a given crystal structure is necessary to evaluate if under 
small deformations of the lattice, the atoms return to their assumed equilibrium 
position, which identify a truly metastable phase [DEI Eli. 

The aim of this study is to evaluate the dynamical stability of bulk ferromagnetic 
Ru and Os in the bet structure with c/a < 1, using DFT-based hrst principles 
calculations. Thus, we determine the structural properties, magnetic moments, 
and phonon dispersion of Ru and Os for the local minimum in the Bain path for 
c/a < 1 and for uniaxial and epitaxial strain in the range of ±5 %. Imaginary 
phonon frequencies were found along several directions of the Brillouin zone for all 
the studied structures. Thus, the results indicate that bulk ferromagnetic Ru and 
Os in the bet structure are dynamically unstable. 

2 Computational details 

The DFT calculations were performed within the framework of the plane-waves pseu¬ 
dopotential approximation (PW-PP), as implemented in the Quantum ESPRESSO 
package in- The exchange-correlation functional was treated with the PBE parametriza- 
tion of the generalized gradient approximation m- The core electrons were replaced 
by scalar relativistic ultrasoft-pseudopotentials from the GBRV database [12], and 
the valence wave functions and charge density were expanded in plane waves with 
a kinetic energy cut-off of 40 and 400 Ry, respectively. 

The integration on the irreducible part of the Brillouin zone was performed on 
a uniform grid of 24x24x24 fc-points, with a small cold smearing of 0.01 Ry [T3] . 

The phonon frequencies were computed by means of linear response theory, using 
the density functional perturbation theory [T3] on a grid of 4x4x4 g-points and 
subsequently interpolated to obtain the full phonon dispersion. 


3 Results and discussion 

The bet crystalline structure belongs to the space group I4/mmm, with an in-plane 
lattice parameter a and an out-of-plane lattice parameter c. For the lattice vectors 
and the hrst Brillouin zone with its high symmetry points, see the appendix A of 
reference [15]. The lattice parameters for Ru and Os in the bet structure were 
obtained from the local minimum in their corresponding Bain paths for c/a < 1. 
In the epitaxial (uniaxial) Bain path, for each a(c/a) the corresponding c/a(a) is 
obtained by a direct minimization of the electronic-total energy, thus making the 
out-of-plane (in-plane) stress disappear [16]. The lattice parameters a and c in the 
local minimum for c/a < 1 will be called abet and Cbct, respectively. 
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Table 1: Lattice parameters abet and (c/ a)bct, energy of the bet structure with respect 
to the hep ground-state {Ebct — Ehep), and magnetic moments (MM) for Ru and Os 
in the bet structure. Previously reported values are included for comparison. 

Metal 

(A) 

(c/(z) j 

^bet E/ficp 

mRy/atom) 

MM 

(/xs/atom) 

Method 

Ref. 

3.24 

0.83 

- 

- 

Exp.'" 

m 

3.25 

0.83 

40 

0.40 

FLAPW 

m 

Ru 3.18 

0.85 

44 

0.62 

FPLO 

0 

3.25 

0.84 

35 

0.60 

FLAPW 

isi 

3.24 

0.84 

38 

0.62 

PW-PP 

This work 


^ 3.30 

3.30 

0.84 

53 

0.50 

FPLO 

0 

0.83 

56 

0.47 

PW-PP 

This work 


^ Data for Ru thin-film prepared by ion beam sputtering on Mo(llO) single¬ 
crystal wafer. 


The calculated structural properties and the magnetic moments of Ru and Os 
in the bet structure are summarized in Table [1] For the lattice parameters, we 
hnd a good agreement with previously reported values [311115116]. Considering the 
calculated energy of the bet structure with respect to the hep ground-state energy 
{Eficp), we determine that Ru and Os have the highest energy difference E^ct — Ehep, 
among the transition metals with an hep ground-state na. The calculated value of 
the magnetic moment for Ru in the bet structure is 0.62 /r^/atom, which is in good 
agreement with recent reports HE]; however, this value is larger than the computed 
value by Watanabe et al. [4] of 0.40 /i^/atom. The obtained magnetic moment for 
Os (0.47 /is/atom) is in good agreement with the value of 0.50 /is/atom reported 
by Schonecker et al. [5]. 

In Fig. [T]we present the epitaxial and uniaxial Bain paths for Ru and Os around 
the local minimum for c/a < 1, including the evolution of the magnetic moment 
for each path. As it can be observed, the epitaxial deformation produces energy 
and magnetic moment changes that are larger than the uniaxial deformation. This 
behavior can be explained considering that the volume per atom is given hy V = 
a^c/2. Thus, for a given material in the bet structure, the deformations on the a 
lattice parameter produce a larger effect on V than deformations in the c lattice 
parameter, and consequently on the total energy and the magnetic moment [19] . It 
is interesting to note that the magnetic moment for Os in the bet structure is more 
sensitive to the epitaxial and uniaxial strains than for Ru in the bet structure. In 
fact, we hnd that the magnetic moment in Os changes from 0 to nearly 0.6 /i^/atom 
in the epitaxial strain range of —4 to -1-4%. This type of behavior is quite interesting 
because of the possibility of developing a strain engineering of the magnetic moment 
in the Os-bet structure. 

To evaluate the dynamical stability of the ferromagnetic Ru- and Os-bet struc- 
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Energy (mRy/atom) 


Epitaxial Bain path Unixial Bain path 



3.06 3.12 3.18 3.24 3.30 3.36 3.42 3.48 0.79 0.81 0.83 0.85 0.87 0.89 


Lattice constant a (A) c/a 

Figure 1: Epitaxial and uniaxial Bain paths for Ru and Os in the bet structure 
around the local minimum for c/a <1. The energy scale is respect to Et^ct- 
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Magnetic moment (pB/atom) 









Figure 2: Phonon dispersion for Ru and Os in the bet structure along the epitaxial 
Bain path plotted in Fig. [H 
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Figure 3: Phonon dispersion for Ru and Os in the bet structure along the uniaxial 
Bain path plotted in Fig. [H 
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tures around the local minimum for c/a < 1, we have calculated the phonon dis¬ 
persion in the epitaxial and uniaxial Bain paths. In Figs. [2] and El we present the 
phonon dispersion for Ru and Os in the bet structure at the local minimum and 
±5% of the epitaxial and uniaxial strain, respectively. In all the cases, we hnd 
imaginary phonon frequencies that clearly indicate the dynamic instability of bulk 
Ru and Os, in the bet structure. It should be noted that the imaginary phonon 
branch in the F ^ X direction reveals an elastic instability related to a negative 
value of the elastic constant C = (cn — Ci 2)/2 [1] |7l HB]. Additionally, the effect of 
the epitaxial and uniaxial strains on the phonon dispersion was evaluated. It was 
determined that neither the epitaxial strain nor the uniaxial strain contribute to the 
stability of the bet structure in both metals, Ru and Os. Furthermore, due to the 
high values of the imaginary phonon frequencies, it is unlikely that the temperature 
effects can reverse that instability. Consequently, ferromagnetic Ru and Os in the 
body centered tetragonal structure with c/a <1 are not truly metastable phases. 

In order to analyze the contribution of the magnetic ordering to the dynamical 
instability of Ru and Os in the body centered tetragonal structure, in addition to 
the ferromagnetic order (FM), we also consider antiferromagnetic order (AF) and 
the non-magnetic system (NM). The total energy calculations for Ru and Os in 
the NM and AF state shows a local minimum along the Bain path for c/a < 1 
(not shown here), with lattice parameters very close to the corresponding FM state. 
Nevertheless, for both metals the energy of the local minimum in the NM and AF 
states are higher than the FM state. For Ru(Os), the AF-FM and NM-FM energy 
differences are 2.1(0.7) and 2.4(0.8) mRy/atom, respectively. 

Considering that the structural parameters for the three studied magnetic states 
(NM, AF and FM) are very similar, we have calculated the phonon dispersions 
for the NM and AF states at the structural lattice parameters of the FM state. 
In this way we are considering only the effect of the magnetic ordering on the 
phonon instabilities. In Fig. E] we show the phonon dispersion for Ru and Os in 
the body centered tetragonal structure with the NM and AF order, including the 
previous FM state for reference. Thus, we found that the NM and AF states are 
dynamically instable for both metals. Furthermore, the values of the imaginary 
phonon frequencies for the NM and AF states are larger than the FM state. These 
results show that the ferromagnetism is not the origin of the instability of Ru and 
Os in the body centered tetragonal structure. 

Finally, we want to emphasize that the predicted dynamical instability of fer¬ 
romagnetic Ru and Os in the bet structure applies to bulk. Thus, the implication 
of this study is that there can be no free-standing bet structures with c/a < 1 for 
these metals. Any bet structures with a c/a < 1 observed for Ru and Os, such as 
the thin dims of Ru on Mo(llO) |3], must be stabilized by an external force. In 
the case of ultrathin dims, the material is stabilized by external edects, such as 
surface energy, substrate strain, low dimensionality and dnite size, among others. 
From a theoretical point of view, a correct drst principles prediction on the stability 
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Figure 4: Phonon dispersion for Ru and Os in the bet structure for the NM, FM, 
and AF order. The AF calculations were performed using the simple tetragonal unit 
cell with a two atom basis. 


and the properties of transition metals in ultrathin dims should explicitly include 
the overlayers and the substrate, and cannot be based only on bulk calculations, 
particularly when the bulk structure of the metal-overlayer is dynamically unstable. 


4 Summary 

We have performed first principles calculations of the phonon dispersion for ferro¬ 
magnetic Ru and Os in the bet structure, to evaluate their dynamical stability. The 
phonon dispersion for the local minimum in the Bain path with c/a < 1 and the 
uniaxial and epitaxial strained structures were analyzed. For both systems, we find 
imaginary frequencies along different directions of the Brillouin zone, which indi¬ 
cated dynamical instability. This result shows that ferromagnetic Ru and Os in the 
body centered tetragonal structure with c/a < 1 are not truly metastable phases. 
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